jordan_product <-
function(x,y){
	c(sum(x*y),y[1]*x[-1] + x[1]*y[-1])
}

